A Self-Consistent Model for Positronium Formation from Helium Atoms 
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The differential and total cross sections for electron capture by positrons from helium atoms are 
calculated using a first-order distorted wave theory satisfying the Coulomb boundary conditions. In 
this formalism a parametric potential is used to describe the electron screening in a consistent and 
realistic manner. The present procedure is self consistent because (i) it satisfies the correct bound- 
ary conditions and post-prior symmetry, and (ii) the potential and the electron binding energies 
appearing in the transition amplitude are consistent with the wave functions describing the colli- 
sion system. The results are compared with the other theories and with the available experimental 
measurements. At the considered range of collision energies, the results agree reasonably well with 
recent experiments and theories. 
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I. INTRODUCTION 

Both positrons (e^) and positronium atoms (Ps) have 
important apphcations in many different branches of 
physics, chemistry and other fields [IHl]- This has moti- 
vated numerous studies of collisions in which electrons in 
atomic or molecular targets are captured by positrons 0- 
[26j . In a collision, each of the various modes into which 
the system under study may be fragmented is called a 
channel. For definiteness, consider the positron-helium 
scattering system, on which the present paper is fo- 
cused. In the entrance channel of the e+ + He system 
the positron impacts on the helium atom. The exit chan- 
nel can be one of a number of possible fragmentations, 
such as elastic or excitation scattering (e+ + He and/or 
e"*" -|- He*), positronium formation (Ps + He^), or ion- 
ization (e"'" + e^ + He*). An exit channel is said to be 
open if the corresponding collision is allowed by all known 
conservation laws, such as energy and momentum conser- 
vation; otherwise, it is said to be closed. Since the transi- 
tion amplitude in general is a function of the energy, the 
occurrence probability for each open channel depends on 
the impact energy of the projectile. The transition prob- 
abilities per unit time, per unit target scatterer and per 
unit of the flux of incident particles with respect to the 
target are called the cross section. Experiments measure 
cross sections, while theoretical studies usually attempt 
to compute them. In positron-helium collisions, since the 
charge transfer process amounts to a four-body prob- 
lem, the theoretical investigation of the rearrangement 
involves the full complexities of the quantum mechani- 
cal four-body problem and as such is too difficult to be 
exactly solved. 

Notwithstanding the difficulty, numerous theoretical 
investigations are found in the literature ll7Hl9l, 33 
[38. . ,40. , ,41j , along with experimental studies 127H32 
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Most of the theories were formulated within a single ac- 
tive electron picture. In this formalism, it is important to 
provide a consistent and realistic treatment of the passive 
electron screening effects, in order to satisfy the Coulomb 
boundary conditions and post-prior equivalency and to 
make the wave functions, binding energies, and Coulomb 
phase factors consistent with each other. 

In the context of K-shell electron capture processes by 
fast protons from multi electron atoms. Decker and Eich- 
ler have discussed the deviations of a few of the standard 
formalisms from these constraints in some detail p^. 
They adopted a parameterized potential due to Green, 
SeUin and Zocher [H, the GSZ potential, to con- 
struct a self-consistent screened boundary-corrected first- 
Born-approximation theory for calculation of the K-shell 
electron capture total cross sections in collisions of pro- 
ton with helium and carbon as well as collision of alpha 
particles with lithium atoms. 

In the present paper a distorted-wave boundary- 
corrected first-order Born (DWBIB) formalism accom- 
panied by the GSZ potential is applied to obtain a satis- 
factory description of positronium formation in positron- 
helium atom collisions at intermediate energies. In this 
model, the single-zeta Hartree-Fock wave function and 
its corresponding binding energy, which show very good 
agreement with those for GSZ potential, describe the ini- 
tial bound state of the electron. Consequently the wave 
function and the corresponding energy in the post or 
prior formalisms are the eigenfunction and eigenenergy 
of the same effective screening potential that appears ex- 
plicitly in the corresponding amplitudes. Thus, in ad- 
dition to satisfying both the Coulomb boundary condi- 
tions an post-prior symmetry, the wave functions, bind- 
ing energies and Coulomb phase factors are consistent 
with each other. 

The plan of the paper is as follows. Section 2 out- 
lines the formalism. Section 3 presents and discusses the 
results and compares them with other theories and ex- 
perimental data. The concluding remarks and a sum- 
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mary comprise the last section. Unless otherwise stated, 
atomic units are used throughout. 



II. THEORY 

Consider a bare ion P with mass Mp and charge Zp 
impingent on an atomic target composed of an active 
electron e and a residual target ion T with masses m and 
Mt, respectively. During the collision, the projectile cap- 
tures the electron. According to the time-independent 
boundary-corrected perturbation formalism of electron 
ca ptu re developed by Toshima et al [H] and Belkic et 
al [43] , the first-order corrected-boundary distorted- wave 
Born(DWBlB) amplitude for such a process has the prior 
and post forms 



A 



[prior) 
DWBIB 



and 



A 



{post) 
DWBIB 



^{xfm-vn\x.) (1) 



(x/i(^/-^r)lx.), (2) 



where and are the asymptotic limits of the dis- 
torting potentials Vi and Vf in initial and final channels. 
In position space, the corresponding distorted wavefunc- 
tions Xi and X} are given by the equalities 



Xi(rT, Rt, R) = 0j(rT) exp(iKi • Rt) 

Zp{Z^-l) 



X exp 



\n{w,R - V, • R) 



(3) 



and 



X/(rp,Rp,R) = 0/(rp)exp(iKf • Rp) 
.Z^{Zp~l) 



X exp I 



ln(v/i?-v/.R) , 



(4) 



in which (j>i{rT) and 0/(rp) are the initial and fi- 
nal bound-state electronic wave functions. In Eqs. ([3]) 
and (j4]), the plane wave functions describe the heavy- 
particle motion in the entrance and exit channels, and 
the phase factors ensure that the overall solutions satisfy 
the proper asymptotic boundary conditions for the dis- 
torting potentials. In these equations, rp and ry are the 
electron coordinates relative to the projectile and target 
nucleus, respectively, Rp is the position vector of the 
center of mass of the Pe subsystem relative to T, Rt 
is a similar vector directed from the center of mass of 
the Te subsystem to P, and R is the internuclear coor- 
dinate vector, directed from T to P. The bare projectile 
and the asymptotic target-ion charges are Zp and Z^, 
respectively, and are the projectile velocity in the 
initial channel and resultant bound subsystem velocity 
in the final channel, respectivley, and and K f are the 
wave vectors describing the relative motions in the two 
channels. 

The Coulomb interaction between the active electron 
and the projectile ion is denoted Vpe(rp) = —Zp/rp, 



with the asymptotic form of V^(R) = —Zp/R. We take 
the internuclear potential to be that of a bare nucleus 
screened by the nonactive electrons represented by frozen 
orbitals. For this purpose, the parameterized potential 
obtained by Green, Sellin and Zocher (GSZ) [ll, is 
adequate because it has the correct behavior near the nu- 
cleus and at infinity and because it yields orbital binding 
energies and wave functions in good ag reement with the 
corresponding Hartree-Fock results . Within a single 
active electron model, using the GSZ potential for a neu- 
tral atom, one can write the projectile-target interaction 
in coordinate representation, Vpt(R) as 



VpAR) = ^[1 + {Zt- IMR)] , (5) 

where Zt is the nuclear charge of the multielectron 
atomic target and fl{R) is given by the expression 



n{R) 



1 



1 - - exp (R/d)] ' 



(6) 



Optimized values for the characteristic parameters H and 
d have been given by Szydlik and Green [i^] for Zt < 18. 
Near the nucleus the projectile ion feels the bare nuclear 
charge Zt, while at large distances the potential behaves 
as Zp/R corresponding to the residual singly charged 
ion. 

With the potential ([5]), the post and prior ampli- 
tudes ll]) and ([2]) satisfy the appropriate boundary condi- 
tions in the initial and final channels. It is easy to see that 
the two amplitudes are identical provided that (tt) and 
(j>f{rp) be the exact solutions of the atomic Schrodinger 
equations associated with the potentials VTei^r) and 
Vpe(rp), respectively. The wave functions, binding ener- 
gies and Coulomb phase factors are moreover consistent 
with each other. Hereafter, given the equivalence be- 
tween the prior and post wavefunctions, we will omit the 
superscripts post and prior in Eqs. ^ and ([2]) and write 
Adwbib to denote the amplitudes. 

The Fourier transform of the VpT potential can be 
written in the form 

TT V TT 



VpT(k) = \l-^ + \l-Zp{ZT-l)V{k), (7) 



where 



Vik) 



dRn{R) sin(fci?). 



(8) 



With the given interaction potentials, Vpe(rp) and 
ypT(R), the transition amplitude for electron capture 
comprises three terms: 

AdWB1B=Ai+A2+A3, (9) 

with the following explicit expressions for the partial am- 
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plitudes Ai , A2 and ^3 in coordinate representation: 
Ai^-Zp j drpdRpx}(rp,Rp,R)-^X»(rT,RT,R) 
A2^ + Zp j drpRpx}(rp,Rp,R)-^X»(i-T,RT,R) 
^3 ^Zp{Zt -I) j drpdRpx}(rp,Rp,R)ffl 



X x»(rT,RT,R) 



(10) 



We now consider the special case Zp ~ 1 and the 
transition from an initial Is hydrogen-like state such as 
(j)i{xT) = Nt exp{—(^TXT) to a final Is hydrogen-like 
state of form (j)f{'x.p) = Np e:xjp{~(pxp), where Np and 
Nt are normalization factors. We apply the Scrodinger 
equation for the final bound subsystem, use the Fourier 
transform analysis, and evaluate the resulting integrals 
to obtain the following expression for the first partial 
amplitude, Ai: 



Ai 



IQtt^NpNtZp{ e/) 

92 



1 



(11) 



where ^ij = Mpm/{Mp + m) is the reduced final-state 
mass of the bound subsystem, and J and K are the mo- 
menta transferred to the target ion and the projectile, 
respectively, during the collision. 

Similarly, using the Fourier transform techniques and 
Lewis integral [45, we derive the following closed form 
for the second partial amplitude A2'- 



A2 = Stt'^NpNtZp 



92 



dCpdCr 



(a2-/?)-i/2ln 



a - (a2-/3)i/2 



(12) 



with 



and 



a = {K^ + (IKt + {J' + CtKp 



(i = {K^ + cDiJ^ + Ct)[v} + id + cl)]. 

Finally, we use Fourier analysis to write the third par- 
tial amplitude ^3 in the form 



A3 = 8tt^NpNtZp{Zt - I] 
where 



1 



dk- 



dCpdCr 
V{k) 



Z(Cp,Ct), (13) 



;k-Kp + C|.)(k + J|2 + C|)- 



(14) 

Before numerical calculations become possible, we have 
to analytically simplify the right-hand side of Eq. (IT51) . 



To this end, notice taken that V{k) is a radial function, 
we use Feynman's identity [ioj to evaluate the angular 
integral over k: 



1 /"^ dx f°° /■+°° 

i{Cp,Ct)^-- I — / dm{R) / dk 

Jo J -00 



Jo 7 



sin(fci?) 

P + 27fc + (5' 
(15) 



where 



and 



= v/J2(l - x)2 + i^22.2 „ 2x{l - a;)J • K, 



5 = {J^+ C|)(l -x) + {K^ + Cp)x. 

Next, we apply Cauchy's residue theorem to the last 
integral on the right-hand side of Eq. (jTS]) to see that 



+00 



dk- 



sm{kR) 



= n-^V^.nn(^fi') 



k^ + 2jk + 6 ^/S^ 

(16) 

With this, we have reduced the four-dimensional inte- 
gral I{(pXt) in Eq. (fT4|) to a two dimensional form. 



Jo 7V" ~ 7 "'0 



dm{R)e-^^/^ sm{-/R), 
(17) 

which can be easily computed numerically, along with 
the derivatives on the right-hand side of Eq. p^ . 

For single-electron capture in impact of positron on 
helium atoms Zt = 2.0, = -0.89648, e/ = -0.25, 
fif = 0.5, Cp = 0.5, Ct = 1.6875, Np = CpTCW^, 
Nt = Ct \/ Ct / TT and the corresponding differential and 
total cross sections are given by 



and 



da 2 Vf 



(^Totai = 27r / (7{9)smede, 
Jo 



(18) 



(19) 



respectively. 

To account for the excited final states, we multiply the 
right-hand side of Eq. (IT^ by 1.202, according to the 
Oppenheimer scaling rule [50| . 



III. RESULTS AND CONCLUSIONS 

We now present the computed differential and total 
cross sections for positronium formation in positron col- 
lisions with helium atoms and compare them with other 
theoretical approaches and experimental data. In order 
to evaluate the third partial transition amplitude, A3, we 
rely on the Gaussian quadrature method to compute the 
two-dimensional integraH(Cp, Ct) in Eq. pT|) . 

Figure [T] shows the angular differential cross sections 
for positronium formation at three incident energies, 100, 
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FIG. 1. Angular distribution of the differential cross sections 
for Ps formation from helium atoms. 
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FIG. 2. Total cross sections for Ps formation in collisions of 
positrons with helium atoms, compared with the results of 
the listed theoretical and experimental studies. 



200 and 500 eV. All three curves display a dip due to 
the cancelation between the contributions of the attrac- 
tive and repulsive interactions to the first-order transition 
amplitude. A similar, unphysical dip has been found in a 
first-order Born treatment of electron capture in the col- 
lision of a proton with hydrogen atoms [51|. Figures (Tf a) 
and (b) compare our results with those of the Coulomb- 
Born approximation (CBA) [52] for three impact ener- 
gies: 100, 200 and 500 eV. In our computation, the dark 
angles, i. e., the angles at which cancelation occurs, are 
31°, 28° and 24°, while the CBA displays dark angles of 
26°, 24° and 23°, respectively. In both treatments, the 
dark angle becomes smaller as the incident energy grows. 



Our results differ significantly from the CBA data at 
100 eV incident energy. Below (above) 30° scattering 
angle our differential cross sections are always smaller 
(larger) than the CBA results. As the impact energy in- 
creases, the difference between the our results and the 
CBA results shrinks. Compared with the difference at 
100 eV, the difference at 200 eV is smaller, and the dif- 
ference at 500 eV is negligible. For energies above 500 eV, 
our results are consistent with the CBA results. Conse- 
quently, the total cross sections obtained from the two 
formalisms converge at high energies. Figure [2] compares 
the total positronium-formation cross sections obtained 
from our formalism with the experimental measurements 
of Diana et al , Fromme et al ^2^ , Fornari et al [s^ , 
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Overton et al 31 1 and Murtagh et al '32'] and also with the 
results of the CBA with Zt = 1.6875, 2 and 7^2^^ 
and of the two-center convergent close-coupling (CCC) 
formalism Is']. As Fig. 2(a) shows, in spite of its sim- 
plicity our formalism yields reasonable agreement with 
experiment at energies above 50 eV. Both the present 
approach and the CBA formalism [H^l are based on per- 
turbation expansion, an approach long known to be un- 
reliable at low scattering energies. 

For the first-order Born approximation to be accurate, 
the higher-order terms in the expansion series must be 
small in comparison with the first-order term. At high 
energies, since two-step scattering processes, such as the 
Thomas double-scattering mechanism, become dominant 
at high energies, the first-order approach yields poor ap- 
proximations to the differential cross sections. All con- 
sidered, we see that first-order methods are only reliable 
in a bracket of energies not so low as to render the pertur- 
bative method inapplicable and not so high as to make 
double-scattering mechanisms dominant. 

Both the CBA and our approach satisfy the correct 
boundary conditions. The two methods are nonetheless 
nonequivalent: while the CBA breaks post-prior symme- 
try, our approach is symmetric, so that the initial and fi- 
nal wavefunctions associated with GSZ potential are con- 
sistent with the scattering potential and bound energies. 
Our approach is hence expected to be more appropri- 
ate to describe experimental processes. As the plots in 
Fig. 2(b) show, however, only small differences separate 
the CBA results from ours. 

More accurate positronium-formation cross sections 
are obtained with the CCC formalism, a more elaborate 
method employing a multiconfigurational wavefunction 
treatment to describe the He ground state. The accu- 
racy of the calculated cross sections is controlled by the 
size of the basis and of the set quantum numbers of the 
included states at each center. Figure HKb) shows good 
agreement between the CCC curve and the experimental 
data. 

Within the energy bracket in which it is reliable, the 
solid line depicting our results in Fig. [2Kb) agrees with the 
CCC line. The difference between the two sets theoretical 
data is less than 0.7% around 50 eV, while deviations of 
approximately 17% separate the CBA and CCC curves 
in the same region. At higher energies the results of the 
three formalisms converge and display good agreement 
with the experimental data. 

Figure 13] compares various theoretical approaches with 
experiment. Cross sections calculated by the CBA with 
Zt — 1.6875, 2 and y/— 2ei [52i] . the correct-boundary 
first-order Born (BIB) jSg] formalism, the distorted-wave 
Born (DW) approximation |39|. the continuum distorted 
wave (CDW) approximation |4C)j and the close-coupling 
method (CC) |4l| are compared with the experiments 
of Fornari et al|30t and Overton et al dJ. The CBA, 
CDW, and CC agree well with the measurments in the 
depicted energy range. The agreement between BIB and 
DW at high energies is fair, but at lower energies those 
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FIG. 3. Total cross sections, compared with the results of the 
CBA [H, BIB [H], DW [lil, CDW H and CC ^ proce- 
dures and with the experimental data from Fornari et al [30l ] 
and Overton et al [3111. 
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FIG. 4. Theoretical cross sections, as in Fig. [3] compared 
with the measurements by Murtagh et al [s^ 



theories deviate from the measurements. 

Our predictions also deviate from experiment at low 
energies, but the agreement improves as the energy in- 
creases. Compared with the results of BIB and DW, our 
results are closer to the experimental data at low ener- 
gies. 

Figure [4] compares the above discussed theories with 
the earlier experimental measurements reported by 
Murtagh et al [12. The CDW, CC, and CBA formalisms 
yield results in very good agreement with experiment 
over the entire energy range. Our results are also in good 
agreement at energies above 75 eV, but the discrepancies 
are sizable at lower energies. The BIB and DW results 
are above the data by Murtagh et al [s^ for all energies 
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FIG. 5. Cross sections for positronium formation in helium 
atoms compared with all available experimental data and with 
the results of the DWA [l7| and CCO 19] formalisms. 



lower than 200 eV, the deviation growing considerably as 
the impact energy decreases. Above 200 eV, all theories 
converge to the experimental data. 

In addition to the CCC formalism, a number of theo- 
retical descriptions of positronium formation in positron- 
helium atom collisions have been recently reportde, on 
the basis of the momentum-space coupled-channel optical 
(CCO) method jJi] and second-order distorted-wave ap- 
proximation (DWA) Figure 5 compares our results 
in the energy range of 50-250 eV with the available results 
from the CCO and DWA formalisms, as well as with the 
experimental data reported by different groups (28l - l32| . 
Our first-order distorted wave results are very close to 
the predictions of the second-order DWA, in which the 
cross sections for Ps formation in the Is orbital has been 
added to 1.66 times the Ps cross sections in the n — 2 
orbitals to yield the total Ps formation cross sections. 

Below 80 eV our results understimate the experimen- 
tal cross sections. At small energies all theories in the 
plot become inaccurate, the the CCO cross sections over- 
shooting the data by Murtagh et al [3^. At relatively 
higher energies, above 80 eV our results agree with the 
other theoretical predictions and with the measurements. 

Figure [S] shows the total positronium-formation cross 
sections in our calculation formalism and all the above- 
mentioned theoretical calculations. Comparison with 
measurements reported by Murtagh et al identifies 
the CDW [40], CC [U and CBA (with Z = 2) as the 
most reliable ones. As before, our data show only fair 
agreement with these three theories, at low energies, and 
good agreement at higher energies. Compared with the 
results of Fojon et al [1] and of Bransden et al [39| , our 
results show substantially superior agreement with ex- 
periment low energies. In the considered range of impact 
energies, our results are compatible with those from the 
more complex CCC formalism. 
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FIG. 6. Total cross sections for Ps formation from helium 
atoms compared the results of the listed theories and with 
the experimental data by Murtagh et al [3^ . 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have applied a first-order distorted 
wave with correct boundary conditions to positronium 
formation from helium atoms. A parameterized potential 
in good agreement with Roothaan-Hartree-Fock wave- 
function and the corresponding binding energy described 
the screening effect of the passive electron on the transi- 
tion amplitude. Our theory is self-consistent, since it sat- 
isfies the correct boundary conditions and the post-prior 
symmetry, and since the interaction potentials, wavefunc- 
tions, binding energies and Coulomb phase factors in the 
formalism are consistent with each other and describe the 
collision and electron capture realistically. 

The chief advantage of our approach is its simplicity. 
The figures in this paper comprehensively compare the 
results of this approach with different sets of experimen- 
tal data and with the results of more elaborate theories. 
They hence set benchmarks monitoring the accuracy of 
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the treatment. Compared with the results of the CBA, 
the angular distribution of our differential cross sections 
are in qualitative agreement even at the smallest ener- 
gies: they differ quantitatively at the smallest energies, 
but the agreement improves as the impact energy grows. 

Our calculated total cross sections are in reasonable 
agreement with other theories, such as the CBA, CDW, 
CC, ceo, and DWA, especially at higher energies. 
Also the results agree well with with the experimen- 
tal data of Fornari et al [13], Overton et al [3l|, and 
Murtagh et al [s^], the agreement with the latter au- 
thors being very good at incident energies above 75 eV. 
As evidenced by comparison with experimental and other 
theoretical data, our procedure becomes less accurate as 
the impact energy is reduced, significant deviations be- 
coming apparent around 50 eV. 

Our relatively simple approach, based on the first- 



order distorted-wave formalism reproduces the results 
of more demanding procedures, based on the second- 
order distorted wave approximation. It can be applied 
to positronium formation in multi-electron atoms, such 
a C, Ne, Na, and Ar, and also to electron caputre from 
atomic shells of multi-electron atoms by bare ion projec- 
tiles at moderate energies. 
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